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ABSTRACT 

The expansion of the Universe leads to a rapid drop in the hydrogen Lya effective optical depth of the inter- 
galactic medium (IGM), f e ff oc (1 +z) 3 ' 8 , between redshifts 4 and 3. Measurements of the temperature evolution of 
the IGM and of the He II opacity both suggest that He II reionizes in this redshift range. We use hydrodynamical 
simulations to show that the observed temperature increase associated with He II reionization leads to a relatively 
sudden decrease in f e ff around the reionization epoch of « 10 per cent. We find clear evidence for such a feature 
in the evolution of f e ff determined from a sample of ~ 1 100 quasars obtained from the SDSS. He II reionization 
starts at redshift w 3.4, and lasts for Az w 0.4. The increase in the IGM temperature also explains the widths of 
hydrogen absorption lines as measured in high-resolution spectra. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: formation — intergalactic 
medium — quasars: absorption lines 



1. INTRODUCTION 

Neutral hydrogen in the intergalactic medium (IGM) reso- 
nantly scatters the flux blueward of the hydrogen Lya emis- 
sion line in quasar spectra (Gunn & Peterson 1965; Bahcall & 
Salpeter 1965; Peebles 1993, §23). The fact that not all the 
flux is absorbed implies that the IGM is very highly ionized. 
The neutral fraction is determined from the balance between 
photo-ionizations, produced by the UV-background radiation 
from galaxies and quasars, and recombinations. For hydrogen, 
the recombination rate depends on temperature oc 1/T 7 , and 
so an increase in temperature will lead to a decrease in the level 
of absorption. Recent evidence suggests that He II reionizes 
around a redshift z ~3-3.5, and the associated temperature in- 
crease should have a measurable effect on the mean absorption. 
In this letter, we compute the optical depth evolution of a sim- 
ulation in which He II reionization heats the IGM at redshift 
3.4. The resulting dip in the evolution of the mean absorption 
matches that recently measured by Bernardi et al. (2002) in the 
Sloan Digital Sky Survey (SDSS) data. 

Hydrodynamical simulations have proved to be very success- 
ful in reproducing many properties of the observed absorption 
(Cen et al. 1994; Zhang, Anninos & Norman 1995; Miralda- 
Escude et al. 1996; Miicket et al. 1996; Hernquist et al. 1996; 
Theuns et al. 1998; Machacek et al. 2000; see e.g. Efstathiou, 
Schaye and Theuns (2000) for a recent review). Most of the 
absorption at redshifts 2-5 arises in modestly over- and under 
dense filamentary and sheetlike structures, leading to a forest of 
Lya -absorption lines in the spectra of quasars (Lynds 1971). 
These structures can be described and understood with simple 
physical models in which the typical size of the absorbers is 
determined by the Jeans length of the photo-heated gas (Bi & 
Davidsen 1997; Schaye 2001). Given that these hydrodynam- 
ical simulations reproduce the data in great detail, we are con- 



fident that they will also allow us to investigate the effect of a 
sudden increase in T on the mean absorption. 

Several lines of recent evidence suggest that He II reionizes 
around z ~ 3-3.5, significantly later than hydrogen. Observa- 
tions of the He II Lya -forest detect a sudden increase in the 
mean He II opacity around z ~ 3 (Reimers et al. 1996; Heap et 
al. 2000; Kriss et al. 2001; Smette et al. 2002). Such a rapid 
increase in the flux of He II ionizing photons is expected to oc- 
cur at the final stages of reionization when individual He III 
bubbles around sources percolate (Gnedin 2000). The result- 
ing hardening of the ionizing background could be responsible 
for the observed jump in the relative abundance of CIV/SilV 
(Songaila & Cowie 1996), although more recent data (Kim et 
al. 2002) apparently do not seem to support such a change (see 
also Dave et al. 1998). Finally, entropy injection associated 
with He II reionization will increase the temperature of the IGM 
(Miralda-Escude & Rees 1994), which has the effect of making 
the Lya -absorption lines broader on average. The gas will 
become nearly isothermal when the change in temperature is 
large. Schaye et al. (2000) detected both signatures in high- 
resolution Keck data, using a method based on the finding of 
Schaye et al. (1999; see also Ricotti et al. 2000 and Bryan & 
Machacek 2000) that the variation of the widths of the absorp- 
tion lines as a function of column density, can be used to mea- 
sure the temperature as a function of the density. Theuns et al. 
(2002a) performed a wavelet analysis of the spectra of several 
high resolution QSO spectra, and found that the data required 
a large increase in T by a factor ~ 2 around z ~ 3.3, consis- 
tent with late He II reionization (Miralda-Escude & Rees 1994). 
The temperature decrease after z ~ 3.3 is also consistent with 
this interpretation (Theuns et al. 2002b). Other heating mecha- 
nisms have been discussed in the literature as well, for example 
galactic winds (Cen & Bryan 2001). Although these may have 
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contributed to the entropy of the high redshift gas, it seems un- 
likely they could result in a sudden change of temperature as 
seems required by the data. 

In the next section, we use hydrodynamical simulations to in- 
vestigate the effect of He II reionization on the evolution of the 
mean effective optical depth, f e ft . Section 3 summarizes the 
method used by Bernardi et al. (2002) to measure f e ff(z) from a 
sample of 1061 SDSS quasar spectra, and section 4 compares 
the data to the simulations. 

2. HYDRODYNAMICAL SIMULATIONS 

We will use hydrodynamical simulations to compute the 
change in the effective optical depth f e ff, i.e. the ratio between 
the observed and emitted fluxes in the Lya -forest region of a 
quasar spectrum, 



exp(-r eff ) = ( 



-observed \ 
-^emitted 



(exp(-r)) . 



(1) 



as the temperature increases during He II reionization. The 
Gunn-Peterson optical depth tqp of a uniform IGM at redshift 
Z is (see Peebles 1993, §23) 



tgp(z) = 



aocxtiH 
H(z) 



(2) 



where <jo = (37r<7r/8) 1/2 /Ao> with aj the Thomson cross- 
section, and/ = 0.416and Ao = 1215.67Athe oscillator strength 
and wavelength of the hydrogen Lya -transition. The Hubble 
constant is H(z), x = «hi/«h is the neutral fraction, and «h is 
the hydrogen number density. In photo-ionization equilibrium, 



X SS (1 + XHeII + 2xHeIIl) 
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where x He ii = «Heii/"H and x He m = «Hern/«H are the abun- 
dances of singly and doubly ionized helium, a ss 3.975 x 
10~ 13 (r/10 4 K)"° 7 cm 3 s _1 is the hydrogen recombination coef- 
ficient, and r is the photo-ionization rate. We have assumed the 
gas to be highly ionized, i<1. Inserting numerical values, and 
assuming a helium abundance of Y = 0.24 by mass, 
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where Q,t,h 2 is the baryon density, f2 m the matter density, and 
Ho = 100/zkms~ 1 Mpc~ I is the current Hubble constant. Since 
the IGM is not uniform, the optical depth will fluctuate, and 
the resulting fluctuating Gunn-Peterson absorption produces 
the Lya -forest. This is what is computed in the simulations. 
Because of the temperature dependence of the optical depth, we 
would naively expect f e s to decrease by ss 1 /g 0J as the temper- 
ature increases by a factor g during He II reionization. We will 
show that the decrease is in fact less than that. 

We have performed a simulation of a vac- 
uum energy dominated, cosmologically flat cold 
dark matter model, with cosmological parameters 
(f7 m ,lWAcr 8 ,F)=(0.3,0.019,0.65,0.9,0.24). The simulations 
are performed using a modified version of HYDRA (Couch- 
man, Thomas & Pearce 1995; Theuns et al. 1998), in which 

8 http://pitto. mib.infn.it/~haardt/cosmology.html 



the gas is photo-heated by an imposed uniform UV-background 
which evolves with redshift and reionizes H I and He I at z ~ 6 
and He n at z ~ 3.4. During reionization of both H I and 
He II we artificially increase the photo-heating rates - i.e. the 
mean energy associated with each photo-ionization - to mimic 
the effects of radiative transfer (Abel & Haehnelt 1999). The 
temperature at the mean density, T{z), of the simulation fits 
the data of Schaye et al. (2000) and Theuns et al. (2002a), in 
particular there is a sudden increase in T at z ~ 3.4 by a fac- 
tor g ss 1.8 as a result of He II reionization. Non-equilibrium 
rates for photo-ionization, cooling and heating are computed 
using the fits in Theuns et al. (1998). The particle masses 
for dark matter and gas are 1.1 x 10 7 and 2.0 x 10 6 M Q re- 
spectively, and the simulation box is 12k' 1 (co-moving) Mpc 
on a side. To illustrate the effect of He II reionization on 
f e ff, we use a similar simulation in which He II reionization 
starts slightly later, at z ~ 3.2. Note that these simulations as- 
sume that the radiation field is spatially uniform at all times, 
see Sokasian, Abel & Hernquist (2002) for a He II reioniza- 
tion simulation that includes some radiative transfer effects. 
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FIG. 1 . — Evolution of the temperature at the mean density T (top panel), 
and the deviation of the effective optical depth f c ff from a power-law (bot- 
tom panel) for simulations in which He II reionization starts at z ~ 3.4 (full 
line and filled circles) and z ~ 3.2 (dashed line and open circles). Data points 
with error bars are the measurements from Schaye et al. (2000). He II reion- 
ization increases T by a factor ~ 1.8 (dotted lines). Away from reionization, 
fat oc (1 +z) c , with ( « 3.8. The increase in T during reionization decreases 
the neutral hydrogen fraction and hence also the mean absorption. This causes 
f c ff to fall below the power law evolution, producing a feature in f c ff/(l +zr 
that is characteristic for the occurrence of He II reionization. 

While the simulation was running, we stored the data required 

to compute mock spectra along random sight lines through 
the simulation box, at many thousands of output times. The 
time interval between individual snapshots is of order of the 
light crossing time through the box. This fine sampling of the 
redshift evolution allows us to take into account the detailed 
evolution of the density, temperature and ionizing background. 
We impose the ionizing background Jg+q(u,z) from QSOs and 
galaxies as computed by Haardt & Madau (1996) and updated 
by Haardt & Madau (2002 s ), but reduced in amplitude by a 
factor of j. Using j = 2.75 for all redshifts 2.5 < 4 reproduces 
the observed evolution of f e g quite well. This UV-background 
is different from the one used during the simulation. However, 
because the hydrogen gas is in photoionisation equilibrium, and 
because the hydrogen photo-heating rate is independent of the 
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hydrogen ionizing flux, such a rescaling works accurately (see 
Theuns et al. 1998). 

Using the mock spectra, we determine the evolution of r e ff, 
binned in redshift intervals of order Az = 0.1, and use bootstrap 
re-sampling to estimate the expected variance for the SDSS 
sample, given the number of QSOs available. This scatter could 
be underestimated due to missing large scale power in the sim- 
ulation box, but such a systematic effect is unlikely to affect the 
relative decrement during and away from He II reionization. 

In Fig. 1, we compare the evolution of T and f e s in the two 
simulations, in which He II reionisation starts at redshifts 3.4 
and 3.2, respectively. The temperature T increases suddenly by 
a factor ~ 1.8 following reionization. The effect of this tem- 
perature increase on f e ff is shown in the bottom panel. Away 
from He II reionization, f e ff oc (1 +zp, with ( m 3.8. The bot- 
tom panel of the figure shows the deviation from this scaling, 
T e s/(J +z) , and illustrates that the temperature increase leads 
to a w 10 percent decrease in f e ff. After Az ss 0.4, f e ff recovers 
to the power-law evolution. 
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FIG. 2. — Evolution of f c ff for the simulation with He II reionization at 
redshift 3.4 (circles). The other curves illustrates two reasons why the change 
in effective optical depth is smaller than naively expected. If peculiar veloci- 
ties are neglected in the computation of the spectra (diamonds), then the mean 
absorption is higher because gravitationally induced peculiar velocities tend 
to increase blending, which in turn decreases the opacity. The curve has been 
off-set vertically to match the full simulation curve above z = 3.5. During He II 
reionization, pressure induced peculiar velocities push gas out of the filaments 
into the voids, and this increases the absorption, and consequently the two 
curves start to deviate. This shows that pressure induced peculiar velocities 
partly counteract the decrease in r e ff . The increased electron density follow- 
ing He II reionization, neglected in computing the dotted curve (squares), has 
a similar effect on f c ff . When both peculiar velocities and extra electrons 
are neglected (dashed curve connecting triangles), then the decrease in f e ff is 
~ (1./1.54 7 ) (parallel dashed lines), close to that expected from the temper- 
ature dependence of the recombination coefficient. 

This decrease is less than the 1/1.8° 7 ~ 30 per cent naively 

expected from the temperature dependence of the hydrogen re- 
combination coefficient. Note however, that f en - does not scale 
linearly with the optical depth, since exp(-f e ff) = (exp(-r)}, 
and that some of the absorbing gas is at densities higher than 
the mean, for which the relative change in temperature is less 
(Schaye et al. 2000). In addition, there are two other physical 
effects that lessen the dependence of f e s on T. The extra elec- 
trons liberated during He II reionization tend to increase the 
neutral hydrogen abundance (Eq.|j). For full He II reionization, 
the increase in f e f t is w 7.5 per cent. Another effect is that the 



filaments start to expand due to the extra heating (see Theuns, 
Schaye and Haehnelt, 2000, Fig (4)). These extra peculiar ve- 
locities also partly counteract the decrease in f e ff . Both effects 
are illustrated in Fig. 2. 
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FIG. 3. — Deviation of the effective optical depth from a power-law evolu- 
tion, r e ff/(l +z) 3 ' 8 , for the SDSS data smoothed on 3000km s~'(symbols with 
error bars) and a hydrodynamical simulation of a ACDM model, in which 
He II reionization starts at z = 3.4 (hashed region). The errors in the SDSS 
data are determined from bootstrap re-sampling the individual spectra. The 
hashed region for the simulation delineates the 20 and 80 per centile of boot- 
strap re-sampled mock samples. The temperature increase associated with 
He II reionization causes f e ff to drop below the power-law evolution in the 
simulation. This characterictic dip matches very well the feature detected in 
the SDSS data. 



3. DETERMINATION OF THE Lya EFFECTIVE OPTICAL DEPTH IN 
THE SDSS DATA 

Bernardi et al. (2002) determine f e ff(z) from 1061 moderate 
resolution (A/ A A « 1800) and intermediate signal-to-noise ra- 
tio (typically between 4 and 10) spectra of QSOs with redshifts 
between 2.75 and 4.3, color selected from the SDSS (Fig. 3). 
There are several crucial steps in the determination of the mean 
absorption from these spectra, such as the determination of the 
shape and amplitude of the continuum in each QSO spectrum, 
and assessment of possible biases introduced by the QSO se- 
lection criterion employed by the SDSS. 

Full details of the procedure can be found in Bernardi et al 
(2002). Briefly, the shape of the continuum in the Lya -forest 
region is determined by solving simultaneously for the evolu- 
tion of the effective optical depth and for the shape of the mean 
continuum, while allowing for the possibility that neither are 
well-fit by featureless power-laws. This is possible, because the 
(intrinsic) shape of a quasar's continuum is a function of rest- 
frame wavelength, whereas f e s is measured at given observed 
wavelength. Two different methods were used, a \ 2 technique, 
and an iterative fitting procedure, and they give results in very 
good agreement with each other. Both techniques allow one 
to take into account the presence of weak emission lines in the 
QSO's continuum, and of deviations in the evolution of the op- 
tical depth from a smooth power-law. 

Bernardi et al. present many tests of their methods. For ex- 
ample, they show that the scatter of individual measurements 
from the mean evolution is uncorrected with rest-frame wave- 
length, which is not the case when the weak emission lines are 
not taken into account. The evolution of f en - determined from 
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the SDSS data is consistent with previous measurements, and 
the excellent statistics allow for the detection of a localized fea- 
ture in feff(z) of the order of 10 per cent. The next section com- 
pares the SDSS data to the hydrodynamical simulation. 

4. RESULTS 

The deviation of ? e g from a power-law f e s/(l + z) 3 ' 8 is plot- 
ted in Fig. 3 for the SDSS sample (symbols with error bars) 
and the simulation for which He II reionization starts at z « 3.4 
(hashed region). With respect to the power-law, the measured 
f e ff shows a distinctive dip around redshift z = 3. 14 ±0.03 of 
depth Af eff = -0.09 ± 0.02 and width Az = 0.09 ± 0.02 (see 
Bernardi et al. 2002 for the determination of the errors), con- 
sistent with the signature of He II reionization. The similarity 
between the data and the simulation is striking. Note that the 
parameters for the simulation were tuned to match the temper- 
ature evolution measured by Schaye et al. (2000), and hence 
the dip in f e g is a prediction of the model. The fact that the 
power-law evolution as measured from the SDSS is the same 
as computed from the simulation, suggests that the ionizing 
background Jg+q(v,z) describes the evolution of the photo- 
ionization rate well. 

The amplitude, shape and duration of the deviation of 
7"eff/(l+z) from a straight line is very similar for the SDSS 
data and the hydrodynamical simulation. The onset of the de- 
crease in f e ff is not correctly modeled in the simulation, because 
the UV-background is assumed to be uniform at all times and 
because the rise in the UV-background does not take into ac- 
count that photons are consumed during the reionization pro- 
cess. This is partly counteracted by the still relatively poor 
sampling Az ~ 0. 1 - worse than in the SDSS data - of the sim- 
ulation outputs, but better sampling would show a much more 
rapid decrease in f e ft in the simulations. In reality, the tempera- 
ture distribution will become inhomogeneous, as some regions 
become ionized - and hence hot - before others do. The spa- 
tial scale of these fluctuations is of the order of the size of a 
new He III region in a singly ionized and homogeneous He II 
Universe ( Miralda-Escude, Haehnelt & Rees 2000) 

16* f JV^^qso. A V Mpc (5) 

1+z ^10 5 V> 10 7 yr l 0.02 j ) P ' K ) 



which corresponds to w 1200km s" 1 or Az = 0.016 at redshift 3 
(/V p h is the luminosity of the QSO in ionizing photons, and Jqso 
its age). Simulations that include radiative transfer and can re- 
solve the Lya -forest are required to model this process more 
realistically. 

In principle, it should be possible to detect the onset of 
He II reionization from the occurrence of such ionized, hot 
He III bubbles. (The 'thermal proximity effect'.) Theuns et 
al's (2002b) wavelet method did not detect such temperature 
fluctuations on scales larger than 5000 kms" 1 . However, they 
only attempted to identify individual He III regions, and did not 
quantify whether the presence of such hotter regions can be de- 
tected statistically. Zaldarriaga (2002) attempted to do just that, 
using the spectrum of QSO 1422+231, but again no signal was 
detected. With more data available now, such a detection may 
be possible. 

In summary: the evolution of the effective optical depth f e ff, 
measured in the SDSS quasar sample by Bernardi et al. (2002), 
deviates suddenly by w 10 per cent at redshift z ~ 3.1 from a 
smooth power-law. The mean absorption determined from a 
hydrodynamical simulation of the IGM, in which He II reion- 
ization starts at redshift 3.4, provides an excellent match to this 
feature in f e ff. In this simulation, the decrease in absorption 
is due to the temperature dependence of the hydrogen recom- 
bination coefficient, and the required temperature change due 
to He II reionization, is in good agreement with the tempera- 
ture measurements from Schaye et al. (2000) and Theuns et al. 
(2002). We consider this strong evidence that He II reionization 
has been detected in the SDSS data. The onset of reioniza- 
tion starts around z ~ 3.4 and percolation, corresponding to the 
maximum deviation of f e s from a power-law, occurs at z ~ 3 . 1 . 
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